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Recent experiments have revealed the role of staggered Dzialoshinskii-Moriya interaction in the 
magnetized phase of an antiferromagnetic spin 1/2 two- leg ladder compound under a uniform mag- 
netic field. We derive a low energy effective field theory describing a magnetized two-leg ladder 
with a weak staggered Dzialoshinskii-Moriya interaction. This theory predicts the persistence of 
the spin gap in the magnetized phase, in contrast to standard two-leg ladders, and the presence of 
bound states in the excitation spectrum. Such bound states are observable in Raman scattering 
measurements. These results are then extended to intermediate Dzialoshinskii-Moriya interaction 
using Exact Diagonalizations. 
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Antiferromagnetic systems with spin gaps have been 
the object of great theoretical and experimental interest 
for the last two decades. In particular, ladder systemsi^ 
have been investigated first theoretically as toy mod- 
els for high-temperature superconductivity'^, and then 
experimentally.— Ladder systems with an even number 
of legs are known to possess a spin gap, while ladders 
with an odd number of legs are gapless, a result reminis- 
cent of the one of Haldane for antiferromagnetic spin-S 
chains^ Upon application of a magnetic field, the gap 
of two-leg ladder systems can be closed, leading to the 
formation of a Luttinger liquid statei^ii In such a state, 
there exists a quasi-long range antiferromagnetic order, 
and inter-ladder exchange leads, at low enough tempera- 
ture, to the formation of a Neel state.® Such field induced 
magnetic ordering is the quasi-one dimensional analogue 
of magnon Bose condensation 1^ 

Recent Nuclear Magnetic Resonance (NMR) 
experiments^- have revealed that the compoundii 
Cu2(C5Hi2N2)2Cl4 (Copper diazacycloheptane to be 
abbreviated CuHpCl), originally believed to be a 
simple spin ladder materiali^iiJiiii, possesses staggered 
Dzialoshinskii-Moriy a^^d^ (DM) interactions along the 
rungs. Although recent Inelastic Neutron Scattering 
experiments suggest that this material is closer to a 
three-dimensional coupled dimer system^^ than to a 
real ladder, it was shown theoretically (using a ladder 
description) that including such a spin anisotropy greatly 
improves the low temperature fits of the magnetization 
curv o^®'^^ and of the anomalous behavior of the specific 
heat in magnetic field. Due to the presence of these 
interactions, the spin gap persists in the magnetized 



phase, and a staggered magnetization is present even 
above the temperature for antiferromagnetic orderingfiS 
The outline of the paper is the following. In Sec. [H 
some basic results on two-leg ladders with staggered 
DM interactions are recalled. Then, in Sec. [TTl we 
develop the field theoretical approach and compute 
the Raman scattering spectra in the Fleury-Loudon 
approximation22i2i. We obtain well defined peaks below 
a continuum. The predictions of the field theoretical 
treatment are then compared to Exact Diagonalization 
(ED) results in Sec. IIIII The existence of peaks is 
confirmed, and is shown to persist to intermediate 
couplings. 



I. INTRODUCTION 

The Hamiltonian of the spin 1/2 Heisenberg two-leg 
ladder reads: 

H ~ J\\^n,p ■ Sn + i^p + JlS„^i ■ S„^2 

n,p— 1,2 n 

where S„.p is a spin 1/2 operator, J±,J\\ > and h = 
g^sH. The staggered DM interaction is: 

ifrnf = 5I(-)"D-(S„4xS„,2) (2) 

n 

with D = Dy. 
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The point group of the Hamiltonian Eqs. ((T])- ([2]) is 
C2v'^ It includes the identity, the reflection symmetry 
w.r.t. a rung (or "parity"), the inversion symmetry w.r.t 
the center of a plaquette and the combined operation 
of the two last ones. The Hamiltonian ((T])~([2]) is also 
invariant under translation of two lattice spacing n — > 
n + 2. 



II. BOSONIZATION STUDY 



A. Weak coupling case 

In this section, we derive an effective field theory de- 
scribing the low-energy, long-wavelength excitations of 
the Hamiltonian Eq. (P)- ([2]), valid in the limit Jj_,D <^ 
J||. The bosonized representation of the Hamiltonian ([1]) 
has been derived for J± < Jn in Refs. |1|28^ 29, 30^31]. We 
follow the notations of Ref. [3l|. Using the bosonized rep- 
resentation for the spin operators {p = 1,2): 



We will derive a weak coupling low-energy theory ap- 
plying bosonization techniqucs^^.i^i^ valid for weak rung 
couplings. Then, in Sec. Ill Bl we will consider the oppo- 
site limit of strong rung exchange and derive an effective 
spin chain Hamiltonian. ^"^'^^ By bosonizing the effective 
spin chain Hamiltonian in the limit of weak DM interac- 
tion, we will show that the low energy theory derived for 
weak rung couplings can be extended to strong rung ex- 
change. In Sec. Ill Cl we analyze the effective field theory, 
and show that it supports breather type excitations, in 
analogy with spin chains with staggered DM interactions 
along the chainsj^i^l 
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where a is a lattice cutoff, and [(j)p{x),dxdpi{x')\ = 
iTTSpyS{x — x'), the Hamiltonian of the two- leg lad- 
der is expressed in terms of the fields (f>a — "^^^^ and 
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where we have written d^dp — nllp. It can beshown^a^ (6»<,) ^ Tr/V2. Using we find the bosonized form of 

that for = the Hamiltonian ^ has gapped exci- Hdj^j: 
tations, with a gap A of order J± and that the fields 
have the following expectation values ((/f^) = tt/^/S and 



Hdm ^ D dx 
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In the presence of a nonzero magnetic field the term 
—h J dxdx4>s/T^ has non-trivial effects on the spectrumii^ 
For D = 0, it induces a commensurate-incommensurate 
transition for h = A.— 1^ In the commensurate state 
{h < A), the system remains gapped and non magne- 
tized, whereas in the incommensurate phase {h > A) 
it is gapless and has a uniform magnetization m > 0. 
If a small D ^ is turned on in the incommensurate 
phase, the existence of a nonzero magnetization allows 
the substitution {dx<t>i,2) ^ m in Eq. ([5]) , giving a term 



mD{cos9i — COS02)- Taking into account the fact that, 
in the incommensurate phase, one still has {9a) = 7r/^/2, 
this term reduces to TOZ)(sin6'a/'\/2) sin(0s/'\/2). The 
gap in the antisymmetric modes is given by^i: Ajjm ^ 

J||(Ji/J,|)3&, giving(sin0Jv^) ~ (Ji/J||)i/(i6^"-4) 
For a magnetic field not too strong, Ka ~ 1/2 so that 
(sin(?jV2)^(Ji/J||)i/4. 

The low energy Hamiltonian describing the behav- 
ior of a magnetized two-leg ladder with staggered 
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Dzialoshinskii-Moriya interaction is therefore: 



dx 
2^ 



Ac 



, i?m(sin(0jV2)) . 

dx sm 

27r V2 



(6) 



Such a description is valid for energy scales much lower 
than \h — A|, and is therefore only applicable in a sys- 
tem which is sufficiently magnetized. Let us briefly dis- 
cuss its symmetries. It is obvious that the Hamilto- 
nian ^ is invariant under parity and translation by two 
sites. From ([3]), translation by one lattice site amounts 
to Oi,2 6*1,2 + TT and 0i,2 4>i,2 + 7i'/2 — Trma, and 
interchange of the chains to 6i 02 and (jji <-> 02- 
The combination of these two transformations leaves the 
Hamiltonian invariant. Such combination amounts to 
9s 9s+ttV2, 9a -9a and <j)s 0s +7r/V2(l - 2ma) 
and 4>a —(^a- It is immediately seen that this combined 
transformation leaves the Hamiltonian ^ invariant. The 
Hamiltonian ([6]) is the one of a sine-Gordon model, and 
its spectrum is massive provided > 1/16. Since in 
the magnetized phased Kg — 1/2, the ladder with stag- 
gered DM interaction remains gapped for all h, as found 
in DMRG and ED^Jl The gap is of order A^m ~ 
J||(L>m/J||)8'^-'/(i6^=-i) - J||(mi:)/J||)4/7. It can be seen 
from Eq. ([3]) that this implies a staggered magnetiza- 
tion with: {S^i) = -(5,^,2) - J||(i:'m/J||)i/(i6if.-i) _ 
J||(ml?/J||)^/''. Since the dual field of 0^ is ordered, the 
correlation functions of the non-uniform magnetization 
along z present an exponential decay instead of the power 
law decay obtained for D = 



B. Strong coupling case 

In this section, we consider again the two-leg ladder, 
Eqs. (Ill)- (2]), but this time, we assume a strong couphng 
limit such that J± ^ J, D. In that case, it is reasonable 
to diagonalize first the rung interaction and derive an 
effective spin chain model to describe the behavior of the 
system under fieldJ^i^ We can write: 
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where the pseudospin r,^ = —1/2 when the n— th rung is 
in the singlet state, and = +1/2 when the n— th rung 
is in the triplet state with 1 + 2 = the absence 
of the DM interaction, the Hamiltonian of the effective 
spin 1/2 chain reads^*^: 



H 



J\\ 



^(■''ri^''^fi+l + ■''n '''n+1 + '^n'^n-hl) 



(8) 



Deriving a strong coupling expansion for Eq. ([2]), we find 
that: 



H 



DM 



^/2^^ 
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Note that substituting simply ([7]) in ([2]) yields a result 
which is incorrect by a factor 1/2. The reason for this 
is that the operators do not annihilate the sin- 

glet state contrarily to their low energy representation in 
terms of the pseudospins. As a result of this, some of the 
contributions to Q are missed if we simply replace in ^ 
the spin operators by their strong coupling expansion ([7]) . 

Let us briefly review the symmetries of the strong cou- 
pling model. First, we notice that the model is invariant 
under a translation of two lattice spacing n n + 2. Sec- 
ond, it is also invariant under inversion n — > —n. Third, 
it is also invariant under a reflexion around the middle of 
a bond n 1— n (or a translation of one lattice site) com- 
bined with a rotation of tt around the z axis. If we turn to 
the two-leg ladder, we notice that that this latter symme- 
try was apparently not present in the original model. It is 
in fact the inversion symmetry around the center of a pla- 
quette in disguised form. Under inversion, singlet states 
change sign (as | tii2> - | T2ii> ^ -| T2ii> - | Tii2)), 
but triplet states do not. As a result, operators that 
transform singlet state into triplet one are odd under in- 
version, and operators diagonal on singlet or triplet states 
are even. Thus, in the effective model, inversion becomes 
the transformation t+ ^'''i-n ^'^'^ '''n i i- ^■ 

a reflection around the middle of a bond combined with 
a rotation of tt around the z axis. 

Eq. ^ shows that the ladder under strong field with 
a small staggered Dzialoshinskii-Moriya interaction be- 
comes equivalent to a XXZ spin chain in a staggered 
magnetic field. In the limit D ^ J± , this problem 
can be treated using bosonizatio n^^'^^ , and the ampli- 
tude of the perturbation can even be found exact lyi^i^i^^ 
Similar results are known to hold in the case of spin 
1/2 chains with a longitudinal staggered Dzialoshinskii- 
Moriya interactioni^i^Si^ Our bosonized theory reads: 



H 



dx 
2^ 



uK i7rUf + -{d,^f 



XD 
aV2 



dx cos 9, (10) 



where u and K can be obtained from ([5]) by Bethe Ansatz 
techniques^ provided that D <C J\\. This is a quan- 
tum sine-Gordon (SG) model^i For K > 1/8, its spec- 
trum contains massive solitons^, which means that the 
original Hamiltonian possesses a spin gap. The mag- 
nitude of the gap wiU be Adm oc J|| (L'/J||)4^/(8^-i) 
provided that K > 1/8. We note that this strong 
coupling bosonized Hamiltonian is of the same form as 
the one derived in Sec. Ill Al indicating that the weak 
and the strong coupling regime are continuously con- 
nected. The formal correspondence between the two 
Hamiltonians is obtained by the canonical transforma- 
tion 6* = f - 9s/V2, (j) = -(j)sV2, K = 2Ks. In terms of 
the sine-Gordon model (fTU)) . the symmetries of the lattice 
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Hamiltonian become continuous translational invariance, 
invariance under the transformation 9{x) — *■ 0{—x) and 
(l){x) -H-x) and ^ 9 + n, D ^ ~D. This lat- 
ter symmetry is a Z2 gauge symmetry resuhing from the 
impossibihty of discerning odd and even sites in the con- 
tinuum hmit. A final symmetry of the bosonized Hamil- 
tonian is 9 ^ —6 and (p —(j). The latter symmetry 
is a consequence of the linear dispersion of excitations 
for 13 = in the scaling limit and is absent in the lat- 
tice system. From the semiclassical treatment^, it is 
seen that parity transforms solitons into antisolitons and 
vice- versa while leaving the breathers invariants. Indeed, 
classically, the soliton (resp. antisoliton) of the Sine- 
Gordon model is given by 9{x) = 4 arctan (exp (2:/^)) 
(resp. 9{x) — 4 arctan (exp (— x/^))) so that the two are 
exchanged by a parity transformation. States odd or even 
under parity can be constructed by linear combination of 
soliton and antisoliton states. Lastly, we note that the 
bound states corresponding to the SG breathers should 
be even under parity as breathers contain an equal num- 
ber of solitons and antisolitons^. 

An interesting aspect of the Hamiltonian (|10p is that 
it can possess breathers (bound states of solitons) as ex- 
cited states provided K > l/4,^^ii^ It can be seen 
from Ref. d that this condition is always satisfied. If 
Mq (X Adm is the mass of a soliton, the masses of the 
breathers will be given by: 

/ TT Tl \ 

Mn = 2Mo sin , (11) 

with n integer and 1 < n < 8A' — 1. Under the trans- 
formation (/i — > — and 9^—9 breathers with an odd n 
are odd, while breathers with an even n are even. 

In the case of a spin chain with a longitudinal 
staggered Dzialoshinskii-Moriya interaction, the effec- 
tive theory is also a sine-Gordon model as Eq. PU)) . 



but with K = 1/2 for not too large applied mag- 
netic field. This leads to a smaller number of breathers 
than in the ladder case^^i^i^'^ Such breathers have 
been be detected in Electron Spin Resonance (ESR) 
measurements^^ii^iiai^i^iSi^i^i^i^ In the case of 
the ladder system however, ESR measurements cannot 
be straightforwardly interpreted. Indeed, ESR absorp- 
tion intensity is proportional to an autocorrelation func- 
tion of {Sf + S2), but this operator is vanishing in the 
strong coupling limit according to ([7]) and in the weak 
coupling limit it is mixing the symmetric and the anti- 
symmetric modes, thus complicating the interpretation of 
the spectra. We suggest in Sec. Ill CI that Raman scatter- 
ing intensity measurements give a direct access to these 
breather modes in the ladder systems. 

Let us now discuss the variation of the number of 
breathers with the magnetization in the case of the lad- 
der. For a magnetization near zero, we have AT < 1 so we 
expect to find 6 breather excitations with masses given 
by (fTTj) . As the magnetization increases, K decreases, 
and for K < 7/8 there will be only 5 breathers. When 
AT = 3/4, i.e. at half the saturation magnetization, there 
are only 4 breathers. As the magnetization increases be- 
yond this value, AT becomes an increasing function of 
the magnetization, and the number of breathers becomes 
again increasing. 

When the magnetization is at half the saturation value, 
a more quantitative analysis of the gap becomes possible. 
Indeed, in that case (r^) = and the exponent AT and 
the velocity u are given as)^*^ 

K ^ 3/4, (12) 

3\/3 ^ 
u = —J\\a. 

Moreover, the amplitude A is given by the integrah ^^'^^i^^ 
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sinh(2^) 



sinhtcosh [(1 - 2^^) i] 2A: 



(13) 



which can be evaluated using Ref.'GO', Eq. (6.3.22) giving: 

;.2 9 /r(2/3)y 

^ ~4^V3 l,r(l/3)j ■ ^^^^ 

This expression yields A ~ 0.51769 in agreement with 
Table I in Ref. [SE Finally, using the expression of the 
gap of the sine Gordon model from Ref. [61I, we obtain 
the quantitative expression of the spin gap: 

3/5 
(15) 



A 



DM 



3\/3r(i/io) / ^^/^r(5/6)r(2/3) 

j|i ~4V^r(3/5) I r(i/6)r(i/3) 




i.e. Adm/J\\ ^ 1.67(_D/J||)3/^. Inserting the result ([T! 
in Eq. (fTTj) leads to the prediction of four breathers with 
energies A„ = 2Adm sin(n7r/10) {n = 1, ... ,4). We no- 
tice that the energy of the lowest breather is smaller than 
the energy of the soliton. If we assume that the spin gap 
coincides with the energy of the lowest breather mode, 
we will find a spin gap Ei = 1.032J|| (D/J||)3/5. This 
formula is in rough agreement with the numerical data 
of Sec. mil To obtain quantitative estimates of the spin 
gap away from half the saturation magnetization, where 
(fT5)) is not applicable anymore, one can use the numerical 
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estimations of the amphtude A obtained in Ref. [H. This 
wiU not be attempted here. 



C. Raman scattering intensity 

In that section, we determine the Raman scattering 
intensity of the two-leg ladder under a magnetic field in 
the limit D, J|| ^ J±. We will assume that the Fleury- 
Loudon approximation2£i^ is applicable, i.e. that the 
frequency of the radiation is much smaller than the Mott 
gap of the materiali^ In that case, the (effective) Raman 
operator acts in the (restricted) Hilbert space of the spin 
configurations and can be written as: 

Or = ^Jijiei ■ nij)(eo ■ ?^^i)S^ • (16) 

where the are unit vectors along the bond directions, 
6/ the polarization vector of incoming light, and eo the 
polarization vector of outgoing light. In principle, the 
strength of the Raman scattering Jij could depend on 
the geometry^ but, as in most studies^, we assume it 
to be constant. 

For simplicity, here we restrict ourselves to simple geo- 
metrical setups, e.g. to the XX, YY and X'Y' geometries 
corresponding to e/ and eo , both along the chains direc- 
tion, both perpendicular to the chains and at 45-degrees 
from the crystal axes (X'=X-|-Y, Y'=X-Y), respectively. 
The corresponding Raman operators can then be written 
as, 

Or^ = 7^(Sn,l • S„+i4 -I- S„^2 ■ S„+i^2) , (17) 
n 

OV = 7 5IS„,i-S„,2, (18) 

n 

or' = ^(of^-oD. (19) 

In the absence of DM interaction, it is expected that 
there is no Raman intensity for frequencies smaller than 
twice the spin gapi ^^i^^'^^i^^ For frequencies equal to ap- 
proximately twice the spin gap, a threshold is obtained 
when interaction between magnons are neglected^i^ 
with the Raman intensity behaving as Ir{lo) (lo — 
2A)^/^. When interactions are added, the thresh- 
old becomes a peak?^ Turning to the ladder with 
Dzialoshinskii-Moriya interaction, in the weak coupling 
case, the threshold at twice the spin gap will be obtained 
from the contribution of antisymmetric modes 9a/4'a to 
the Raman intensityj^i The symmetric modes will give 
extra peaks and thresholds but at energy scales small 
compared to the gap. In the following, we will compute 
the Raman intensity for frequencies smaller than 2J±. In 
that case, we can use the strong coupling theory for the 
calculation. We will consider two cases, that of a field 
parallel to the legs in Sec. Ill CTl and then that of a field 
parallel to the rungs in Sec. Ill C 21 



1. Electric field along the legs 

In the strong coupling limit, we can use Eq. ([7]) to 
rewrite the Raman operator pT]) as, 

n 

+ « + l/2)(r„V + 1/2)), (20) 
yielding the expression: 

Op. =7 l-const(^l) 

•^11 

where H is the full strong coupling effective spin chain 
Hamiltonian. In the computation of the Raman scat- 
tering intensity^, we had already noticed that the time 
independent terms would drop out from the calculation. 
Therefore, we can consider the effective Raman operator: 

Of"" - l{ho Y.<- Hdm)/J\\ , (22) 

n 

where ho — J±~ J\\/'2- Going to the continuum limit, 
it is convenient to introduce a density for the Raman op- 
erator, 0^^{x) so that Of^ = JO^^{x)dx. Using the 
density of Raman operator, and the standard definition 
of Raman scattering intensity, we finally find that: 

/^,ll(c.) = ^|(n,P„ = 0|Of^(0)|0)p<S(c.-i?„) (23) 

n 

where |0) is the ground state (GS) of the system and 
\n,Pn) is the n—th excited state with P„ its total mo- 
mentum. The Raman operator being given by a sum over 
all sites, it is easy to show by Fourier transformation that 
Raman scattering is conserving photon momenta. This 
explains why only the states with zero momenta (i.e. hav- 
ing the same momentum as the ground state) contribute 
to the sum. 

To proceed with the calculation, we need the bosonized 
expression of the density of Raman operator. It is easy 
to show that it reads: 

O^^(.) = -^(^cos0(.)-^5.0) (24) 

The Raman intensity can now be found by applying the 
Form Factor metho d^^'^^i^° . The form factors of an oper- 
ator are simply the matrix elements of this operator be- 
tween the ground state and an excited state. In the sine- 
Gordon model, the excited states are formed of solitons, 
antisolitons and breathers with given rapidities. Because 
in one dimensional integrable models collisions between 
particles do not lead to production of new particles^, the 
form factors are derived from the sole knowledge of the 
exact S-matrices of the sine-Gordon modeL ^^'^^i^^ 

In our case, the relevant form factors have been al- 
ready computed by other authors in slightly different 
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contexts™'^'*. Using the symmetries of the operators in 
Eq. pi)) . it is possible to further simphfy the expres- 
sion of the Raman intensity. The sine-Gordon Hamilto- 
nian is invariant under the simultaneous transformations 
6^—9 and cj) ^ —cj), and its eigenstates can be classi- 
fied according to their parity under this transformation. 
We note that the operator cos 9 is even under such trans- 
formation, whereas the operator dx4> is odd. Therefore, 
the non- vanishing matrix elements of cos 9 are between 
the ground state and even states, whereas those of dxtj) 
are between the ground state and the odd states. An 
immediate consequence is that: 



(25) 



where: 



4:1, (-) 



4:1, (-) 



^ \{n\d,m\'Siio-E„) (26) 

^1 \n)odd 

(^77f) E \{n\coa9\0)\'S{iv-E^) 

\ II / \n)even 



The intensity Ij^^^{uj) has been computed previously 

in Ref.^^ in the context of spin chains with staggered 
Dzialoshinskii-Moriya interactions. ^From the results of 
Ref. ,7^, /jj'^ II (w) contains delta-function peaks at the fre- 
quencies of even breathers, and thresholds at frequencies 
equal to twice the soliton frequency, and to frequencies 
equal to the sum of the masses of two breathers of iden- 
tical parity. The lowest threshold is therefore obtained 
at the frequency uj = 2Mi . We notice that this threshold 
is below the energies of the breathers of mass AI3 and 
M4. Therefore, peaks are obtained in the continuum. 
Of course, in a more realistic model, the conservation 
laws that make the sine Gordon model integrable would 
be absent, and one should expect that the peaks would 
be broadened by coupling with the continuum. If the 
deviations of the lattice model from the continuum sine- 
Gordon model are sufficiently important (for instance if 
the gap is so large that the continuum limit is not justi- 
fied), they may even lead to the complete disappearance 
of the two peaks. 

Turning to the intensity /jj°||(ti;), we use the equation 
of motion for the field 9, dt9 — ^dx4>, to relate it to 
the electrical conductivity of a one dimensional Mott 
insulator Ref. |75i, it was shown that if the sine- 
Gordon model describing the low energy charge excita- 
tions of the Mott insulator possesses breathers in its spec- 
trum, the conductivity of the Mott insulator has delta 
peaks at the breather frequencies. In the context of the 
Mott insulator, these peaks were interpreted as exciton 
lines."S Translating the results of Ref. [zE we conclude 
immediately that (lu) possesses delta function peaks 
at the frequencies of the odd breathers. Moreover, there 
are thresholds at frequencies equal to twice the soliton 



mass and to the sum of two breather masses of differ- 
ent parities. The first threshold in 1^°^^^ (to) is thus at the 
frequency lu — Mi + AI2 ■ 

Combining the two contributions, both the odd and 
the even breathers yield peaks in the Raman intensity. 
Moreover, the Raman intensity exhibits thresholds at all 
frequencies equal to the sum of two breather masses or 
to twice the soliton mass. A qualitative sketch of the 
Raman intensity is visible on the figure [TJ 




(a/2M, 







FIG. 1: A sketch of the Raman intensity for polarization par- 
allel to the legs. The arrows represent the position of the the 
breathers (5-peaks. The intensity in the peak, proportional 
to its height is decreasing with the breather index. The con- 
tinuum starts at a frequency equal to twice the mass of the 
lightest breather. Above the threshold, resonances are ex- 
pected instead of J- functions (beyond a SG analysis). 

In the special case of /leS = 0, the intensity /^"^ is 
vanishing. Only the even breathers will contribute to the 
Raman spectrum. 

The weight of the delta peaks can be computed by 
the Form factor expansion. ^^'^^'^^ From the general form 
factor expansion, we have that: 



iR{Lo) = Y,mdR\B^m\ 



, 5{uj - En) 

En 



(27) 



where |i?„(0)) is the state with one i?„ breather of 
momentum zero. En the energy of the breather and 
{Q\Oii\Bn{Q)) is the one breather form factor. Note that 
this formula a priori applies to both || or _L setups cor- 
responding to the XX and YY (see below) polarizations. 

Our task is thus to compute the one breather form 
factors for the two operators dt9 and cos6'. Let us first 
sketch the idea behind the calculation. The axioms of 
Form- Factor theory state that the existence of breathers 
result in poles in the n-particle form-factor (considered as 
complex function of the rapidities of the particles) when 
the difference of two rapidities is equal to a (purely imag- 
inary) fusion angleiSi^S The residue at the pole is pro- 
portional to the {n — 2) particle, one bound state form- 
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factor. Therefore, the Form-factor for a single sohton 
bound-state wiU be obtained from the residue of the two 
particle form factor at the appropriate fusion angle. 

To be more precise, in the case of the sine-Gordon 
model, one has for the soliton-antisoliton form factor of 
an operator O: 



where (p = (pi — (p2, And the one of e'^ by: 



. ,g. , cosh((^/2)e^(*— 

\e"'\ipi,ip2) = . ' . ^ F,nin{(p), (32) 



sinh X{iTT — ip) 



R 



lU 



where: 



(28) 



Where: 



^/2 



|Res5'(6i)| 



1/2 



(29) 



with = (0|O|i3) the form factor of the breather, u the cor- 
responding fusion angle, and S the S'— matrix for soliton- 
antisoliton collision in the sine-Gordon model. In the 
case of the sine-Gordon model^I, the fc— th fusion angle 
is given byu = 7r(l — A:/A), with A = 8K— 1. The residue 
of the S'-matrix is given by: 



exp 



sm 



dt sinh 1(1 - 
t cosh I sinh 



^(^-»y^) 

sinht 



(33) 



|ResS(0)|,^„(,_,/^) = 2cot ( ^ ) [] cot 



Trfc 



The soliton-antisoliton form-factor for dtO is given by: 



{Q\dte\ipi,ip2) = —cos\ ^ 



cosh ^(iTr — tf) 

cosh ^(iTT — Lp) 



Fymniv), (31) 



As a function of the fusion angle, the energy of the 
breather is 2Mcos(w/2) = 2Af sin(fc7r/2A). The Form 
factor of dtO has poles only for odd n whereas the Form 
factor of e** has poles for both odd and even k. Hov- 
ever, the poles at odd k are the contribution of svaO and 
therefore, only the poles for even k must be considered 
in order to obtain the contribution of the cos 6 term to 
the Raman intensity. We find: 



|(0|at(?|B2„+i(0))|2 =Cisin 



7r(2n+ 1) 
A 



2n 



Jltan^lg'jexp 



-2 



dt sinh 1(1 -i) sinh 



t cosh I sinh (^) 



2 / (2K+l)f 
2A 



sinh t 



, (34) 



recovering the result in Ref. |75„ and: 



|(O|cos0|B2„(O))|^ = Casin" 



9 / 27rn\ /"Kn\ 

in"^ tan 



/Trnx'fr' , /-^^ d< sinh |(l - i) sinh^ (|^) 

tan — II tan^ — exp -2 / ^ , 

V A y ^J- \2Xj ^[ Jo t cosh I sinh (^) sinht J' 



(35) 



where Ci and C2 are independent of n. In the case of 
if = 3/4, we have 4 breathers. The weight of the peak as- 
sociated with the first breather is obtained numerically as 
O.345C1. For the peak associated with the third breather 
the weight is only O.OII6C1 i.e. about 3% of the weight 
of the first breather. Similarly, the weight of the second 
breather, we is 0. 03676*2 and the weight of the fourth 
breather is O.OO58C2 i. e. about 16% of the intensity 
associated of the second breather. Note that the ratio 
C2 / Ci depends not only on the Luttinger parameter but 
also on the ratio -D/Jy. 



2. Electric field along the rungs 

In the case of an electric field parallel to the rungs, 
the Raman operator (fT8|) can be rewritten, in the strong 
coupling limit, as: 

Or=7E(-"-i) ■ (36) 

In the absence of the Dzialoshinskii-Moriya interaction, 
this term is proportional to the total magnetization, 
which leads to a vanishing Raman intensity for w 7^ 0. 
However, this result is only valid for frequencies lu <ti Ji_ 
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since the strong coupling approximation eliminates the 
states = 0, 1 that are placed at energies '-^ Jj_ above 
the triplet and the singlet. From the results of Ref. IbsI . 
we know that for lu equal to twice the spin gap, a pho- 
ton can create a pair of magnons, resulting in a non-zero 
Raman intensity. Therefore, in a two-leg ladder with 
D = the Raman intensity should be zero in an interval 
[0,2(J^- Jjl)]. 

When we turn on D, the total magnetization is not 
anymore a good quantum number, and we expect a non- 
zero Raman intensity also for energies small compared 
to J±. Using bosonization, we can express the Raman 
intensity as: 

IrA^) ^^J2^n\mn)\'SiP^)Si^-E,,) (37) 

n 

It is obvious that Ir^± (lo) oc /j^ y {lo) and therefore 
only contains contributions from odd states. Moreover, 
the continuum starts at the frequency Mi + M2 which is 
larger than in the case of a polarization along the legs. 
The behavior of In^^i'-o) is sketched on Fig. [21 Again, as 
we noted in Sec. Ill C 1[ the breather B-^ has a higher en- 
ergy than the lowest two-Si breather excitation. So the 
breather may be unstable in the more realistic case, 
and may appear as a broad resonance or be completely 
absent. 




0.5 1 1.5 2 



(B/2M0 

FIG. 2: A sketch of the Raman intensity for polarization par- 
allel to the rungs. The arrows represent the position of the the 
breathers (5-peaks. The intensity in the peak, proportional to 
its height is decreasing with the breather index. In contrast 
to the case of polarization along the legs, only odd breathers 
contribute to the Raman intensity. The continuum starts at 
a frequency equal to the sum of the masses of the lightest 
odd and even breathers, which is higher than in the case of 
polarization along the legs. 

In the case of X'Y' polarization, it is obvious that the 
resulting Raman operator being a linear combination of 
and O)^^— , the resulting intensity will be qualita- 
tively similar to /fl' y shown on Fig. [T] 



III. NUMERICAL EXACT DIAGONALIZATION 
RESULTS 

A. Introduction 

Since, strictly speaking, the previous analytical ap- 
proaches are only valid in the limit J|| <C Jl, it 
is of great interest to benefit also from a complemen- 
tary numerical approach valid within a broader param- 
eter range. In this Section, coming back to the original 
ladder model of Eq. (IT])-(l2]) (again with alternating DM 
vector and external magnetic field perpendicular to it), 
we perform extensive Exact Diagonalizations (ED) on fi- 
nite periodic 2 x L clusters (the length L ranging from 4 
to 14) using Lanczos (for GS static and dynamical prop- 
erties) or Davidson (for computing low-energy spectra) 
algorithms. 

To illustrate our results, we use a typical parameter 
J|[ = 0.2 J_L consistent with a (crude) theoretical descrip- 
tion of the molecular solid CuHpCl^'^ . When not speci- 
fied otherwise, the largest parameter Jj_ sets the energy 
scale. Fortunately, for such a small ratio, finite 

size effects remain small and results are quite reliable, 
in particular for D > 0.1. Note that finite size scal- 
ing can also often be performed to improve further the 
accuracy of the computation. Interestingly, accuracy is 
getting better for increasing D (at finite fields) i.e. pre- 
cisely in the range of parameters for which the analytic 
treatment is becoming unreliable. The bosonization and 
ED techniques are therefore complementary. Note also 
that, despite the lack of numerical accuracy in that limit, 
we have found that the ED results are consistent with the 
bosonization when D ^ 0, as discussed later. 

As mentioned in previous studiesi^i^^, the DM term 
has been shown to produce a smooth magnetization 
curve, even for finite clusters, and a gap above hd- For 
convenience we reproduce here in Fig. and (b) the 
(reduced) magnetization curve and the gap, respectively, 
as a function of magnetic field for D = 0.1 and D = 0.15, 
values of the DM interaction that we shall consider later 
on, and also for D — 0, for comparison. 

In the following, for simplicity, we shall focus on the 
special field such that m = O.bmsat- From direct inspec- 
tion of the magnetization curve, one gets h ~ I.1J_l— , 
almost independently of the value of D since all curves 
seem to cross at his point. Note that, for such a field, 
bosonization predicts K = 3/4 and four breathers. 

B. Low energy spectra: solitons and bound states 

Prior to the calculation of the Raman spectrum, it is 
interesting and useful to first explore the structure of 
the low energy spectrum. The clusters being periodic 
(ring geometry) , one can take full advantage of the lattice 
translation symmetry (the unit cell contains two rungs) 
so that eigenstates are labelled according to their mo- 
mentum along the legs. We shall focus here on the q = 
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FIG. 3: (Color online) Magnetization (a) and gap (b) vs (re- 
duced) magnetic field for finite D = 0.1 and D — 0.15, com- 
pared to the D = case, computed on a 2 x 12 ladder. Note 
that the h — hc,i singularity at _D = disappears for finite D. 
(The behavior of the gap at _D = is schematic) . All energies 
are in units of Jx. 



sector which is relevant in an optical experiment. As a 
consequence of the C2v symmetry, the g = states can 
be split in four separate symmetry sectors (p, i) = (±, ±), 
corresponding to ±1 characters, respectively. Although 
the Raman spectrum will only involve finite excitations 
belonging to the (-f , -f ) sector (the Raman operator bears 
the full lattice symmetry) , it is useful to consider the four 
possible symmetry sectors. Before presenting our results 
about the four spectra, it is necessary to briefly discuss 
the possible connection of the above quantum numbers 
p and i defined for the lattice model with symmetries of 
the bosonized effective model. First, we note that parity 
corresponds also to a symmetry of the bosonized hamil- 
tonian and of the Sine-Gordon model. In contrast, un- 
der the above inversion symmetry, the bosonized Hamil- 
tonian is not invariant by rather transforms according 
to _D — > —D and 6 ^ 9 + tt corresponding to a dis- 
crete gauge symmetry (the notion of "even" and "odd" 
sites disappears in the long wavelength limit). There- 
fore, the excitations of the bosonized Hamiltonian can 
only be classified according to p — ±. As we discussed in 
Sec. Ill B1 the states in the field theory can only be charac- 
terized by their parity. We recall that we have found that 
linear combinations of one soliton and one antisoliton 
states could be either odd or even under parity, whereas 



breathers had to be even under parity. These simple 
symmetry-based considerations will in fact be quite use- 
ful for a precise identification of the various eigenstates 
found in the numerics. 
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FIG. 4: (Color online) Low-energy spectra in the four differ- 
ent (p, i) symmetry sectors for 2 x L ladders, L ranging from 
L = 4 to L = 14 according to color/symbol codes shown on 
the plots, and for D = 0.1 (a) and D = 0.15 (b). All ener- 
gies are measured from the GS energy {E — 0). Full (LA- 
PACK) diagonalization (small size L = 4, 6), Lanczos and/or 
Davidson ED {L = 8, 10, 12, 14) have been performed*^. The 
tentative onsets of the 2-particle continua are indicated by 
arrows. 



We have computed the low-energy spectrum in each 
symmetry sector of 2 x L periodic ladders with L — 
4, 6, 8, 10, 12 and 14. Results for D = 0.1 and D = 0.15 
are shown on Fig. [D At low energy, well-defined dis- 
crete energy levels clearly saturate to constant values for 
increasing system sizes. Remarkably, we observe that 
these states remain separated from the rest of the spec- 
trum which can be characterized as a "continuum" from 
the accumulation of levels appearing with increasing clus- 
ter size. The discrete levels then correspond either to the 
soliton or to the breathers of the previous Sections to 
which we now try to make a precise assignment. To do 
so, we first notice that two of these levels in the (-I-, -I-) 
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and (— , +) sectors are almost degenerate and might cor- 
respond to the even and odd parity sohton excitations 
(see above). As shown in Fig. El^a), a finite size scahng 
analysis indeed proves that the ratio of these frequencies 
converge to 1 in the thermodynamic limit so that these 
levels can indeed be assigned to the solitons in Fig. IH 
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FIG. 5: (Color online) Ratio of the two quasi-degenerate ex- 
cited states in the (-I-, -I-) and ( — ,+) sectors (labelled as "S" 
in Fig. m as a function of inverse system size for different 
values of D (as indicated on the plot). 

After having identified the soliton level (of energy Mo), 
we can also safely identify the lowest energy excitation 
(which is well separated from the rest of the spectrum) 
as the Bi breather of the SG model. We now tentatively 
assign the remaining discrete levels as B2 , and B^ , as 
shown in Fig. 31 and provide a refined analysis that vali- 
dates this initial guess. First, we note that all suspected 
breathers indeed correspond to p = + as expected from 
the SG theory. Secondly, we use the soliton and Bi lev- 
els as "references" for the identification of the rest of the 
low-energy spectrum by performing the following analy- 
sis: from Eq. pd|) of the previous bosonization approach, 
one expects that the following combinations of the mass 
ratios. 



.2 — 



Ri.: 



( Ml 



( M2 



V2Mo 
IM3 
3Afi 3 V2M0 



V2M1 
4 / A/i 



(38) 



converge to 1 in the limit — > 0. Such relations can then 
be considered as useful "consistency checks" of our initial 
assignments. We have plotted R\,2 and i?i_3 in Fig. [6] as 
a function of D and for different systems sizes. These 
plots are indeed consistent with the fact that R\;i — > 1 
and _Ri^3 ^ 1 for L — > 00 and (then) Z? — > 0, validating 
the assignments of the Bi and breathers. As far as 
B4 is concerned, its identification is more risky but still 
very plausible, at least for D > 0.1. Lastly, we note that, 
as seen in Fig. 01 simple kinematic rules can be used to 
combine two of the above excitations (summing up their 
energies) and lead to the right order of magnitude for 



the two-particle thresholds (shown by arrows in Fig. |4]) 
as well as the correct quantum numbers.. 




0—0 L= 


6 


□ — □ L- 


8 




10 




12 


< — <i L= 


14 


L= 


16 



1 


(b) . 














1 





0.15 



0.05 



0.1 



0.15 ' 



D 



FIG. 6; (Color online) i?i,2 (a) and Ri.s, (b) of Eq. JM]) as 
a function of D for different systems sizes (as indicated on 
the plot). The expected L — > 00 extrapolation is tentatively 
sketched as a dashed line. The magnetic field is such that 
m = O.Smsat and D is measured in units of ,J±. 

As a last check, we have considered another (naive) 
alternative for the Bj, level as the onset of a 2-particle 
continuum made of independent Bi excitations. How- 
ever, the extrapolation of the M3/2M1 ratio (not shown) 
in the thermodynamic limit is inconsistent with M^/2Mi 
being exactly 1 for all values of D. This gives additional 
credit to our previous assignments and suggests the ab- 
sence of a continuum at 2Mi. 



C. Raman scattering and spectral weight 

We now turn to the calculation of the Raman spectrum 
and investigate the effect of the DM interaction. 

It is physically instructive to write the (T = 0) Raman 
spectral function in the Lehmann's representation. 



1 



H + E0 + iO+ 



J2\{p\OBm'K^-iEp-Eo))., 



Ofl,|0)(39) 
(40) 



in terms of the physical excited states \p) of excitation 
energies Ep — Eq. When the Hamiltonian is isotropic 
in spin space, since the Raman operator is SU(2) spin- 
symmetric, it physically describes singlet excitations or, 
more picturially, "double magnon" excitations. In our 
case, because of the Hamiltonian anisotropy, the total 
spin is no longer a good quantum number and the only 
remaining symmetries (apart from the obvious conser- 
vation of the full polarisation m) are the space group 
symmetries discussed above. Operators involving light 
scattering (like the Raman operator) are translational in- 
variant (on the scale of the lattice spacing) and so involve 
only g = excitations. In addition, for light in-coming 
and out-going polarizations along high-symmetry direc- 
tions (see later) , the Raman operator bears the full lattice 
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symmetry so that the Lehmann sum only contains (+, +) 
excited states. However, the actual spectral weight might 
depend strongly on the polarisations. Note that the latter 
selection rule exclude the B2 breather from the Raman 
spectrum. 
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FIG. 7: Raman spectra, computed on a L = 14 cluster, for 
the three different polarizations XX, YY and X'Y' (see text) 
and for D — 0.15 and a magnetic field such that m = O.Smsat • 
The spurious ui = peak has been subtracted (see text). Note 
that a small imaginary part e = 0.01 in the frequency has been 
used to give a small width to the low-energy delta-peaks. The 
corresponding weights of these peaks are shown in Fig. [8] as a 
function of system size. Insets: relative weight of the S state 
with respect to the Bi one vs D, compatible with a vanishing 
S intensity as D — > 



We have used a standard Lanczos technique to com- 
pute Ir{llj) on finite size clusters up to L = 14. Within 
the Lanczos scheme, dynamical correlations can be con- 
veniently expressed as simple continued fractions starting 
from Eq. ([55]) . The spectral weights can also be exactly 
computed from the coefficients of the continued fractions. 
We have checked on small sizes that Lanczos and full di- 
agonalizations do give identical results. We have also 
checked that, in the absence of magnetic field and DM 
{h = and D — 0), our results are compatible with 
Ref. [sj and can be simply interpreted as a large weight 
located at twice the spin gap (the Raman operator is a 
spin singlet and light will create 2 triplet excitations as 
stated above). 

Typical Raman low-energy spectra are shown in Fig. [7] 
for the same parameters as before. Note that the finite 
GS expectation values of the Raman operators produce a 
peak at zero- frequency. To clarify the plots, we have sub- 
tracted this trivial contribution (which experimentally is 
meaningless) and normalized the rest of the spectra so 
that the total integrated weight is unity. Interestingly, 
signatures of the soliton and of the Bi, and even ^4 
breathers (all in the (-I-, -I-) sector, see Fig. 2]) can clearly 



0.01 



(a) XX 



G — e — e- 



o— ^ — ^ 



0.05 



0.15 



0.2 



CUD 



0.01 



0.0001 



0.01 



0.0001 



(b) YY 



o' o 6 6 





FIG. 8: (Color online) Finite-size scaling of the Raman spec- 
tral weights of the soliton and the breathers at D = 0.15 and 
for a magnetic field such that m = 0.5msat. AH light polari- 
sations XX, YY and X'Y' are shown. A logarithmic scale is 
used to reveal the tiny weights in the YY polarisation. For 
XX and X'Y' polarisations a sizeable fraction of the overall 
weight (normalized to 1) is located at higher energies (not 
shown). 



be seen at low frequency in the XX and X'Y' spectra. 
However, for YY polarisation, only the Bi peak is size- 
able. Note also that the low-energy integrated weight 
(let us say, up to w ~ Jl = 1) depends significantly on 
the light polarisation. This is directly linked to a large 
transfer of weight to/from the high-energy region, namely 
above ^ 2Jj_ (not shown), when changing the light po- 
larisation. We have also carried out a finite size scaling 
analysis of the weights associated to each of the discrete 
levels, as shown in Fig.[51 We also find that (for the same 
parameters as before) the relative weight of the S state 
with respect to e.g. the BI breather is 0.35, 0.23 and 0.11 
for D = 0.15, D = 0.1 and D = 0.05 respectively, e.g. 
in the XX polarization (see inset of Fig. [7|). This shows 
that the relative weights of the peaks are rather sensitive 
to the parameter D. Our results show small finite size 
effects and hence strongly suggest that all the sizeable 
peaks in Fig. [7] survive in the thermodynamic limit. 

Lastly, we come back to the issue of the gap which is 
plotted in Fig.[3^b). From the previous analysis this gap 
corresponds in fact to the excitation energy Mi of the 
Bi breather. Its behavior as a function of the DM inter- 
action is shown in Fig. [5] and compared to the bosoniza- 
tion prediction ~ 1.03J||(£'/ Jy)'^/^. Apart from a slightly 
different multiplicative factor, the agreement between 
bosonization and (extrapolated) numerical results is ex- 
cellent. 
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FIG. 9: (Color online) Evolution of the gap (Mi) as a func- 
tion of D for a magnetic field such that m = 0.5 maat (log 
scale). The results for several lengths are reported and the 
extrapolation for infinite size is consistent with a D^^^ be- 
havior (continuous line). The bosonization results shown as 
a dashed line differs only by a multiplicative factor of order 
1.25. All energies are in units of J±. 



cos 9 cos 24> and cos 6 sin 2(f), the Raman operator will con- 
tain soliton creation/annihilation operators. This sce- 
nario could explain the presence of the peak, while pre- 
serving the identification of the state S on Fig. 3] as a 
soliton. Moreover, this scenario is also compatible with 
the decrease of the weight in the peak S as Z) is de- 
creased. Note however that, in this appealing scenario, a 
question remains: why is the B2 breather odd w.r.t. the 
inversion symmetry in contrast to the other breathers ? 
Solving this puzzle unfortunately requires to go beyond 
the sine-Gordon model ^TU\i which does not exhibit such 
a symmetry. 

The third possibility is that the antisymmetric modes 
discussed in Sec. Ill Al are also contributing to the spec- 
trum of Fig. 0] when D ~ Jy even at energy scales small 
compared to J±. In that case, the continuum field the- 
ory method of Sec. Ill Bl would be useless to understand 
the spectra as it implicitely assumes that antisymmet- 
ric modes are excited only for energies comparable with 
J±. The fuller theory of Sec. Ill Al would also be of lim- 
ited usefulness as the resulting generalized sine-Gordon 
model discussed in III Al is non-integrable and thus does 
not permit to locate breather type excitations. 



D. Discussion 

To conclude this Section, we would like to discuss in 
more details the comparison of the Raman spectra ob- 
tained from exact diagonalization and from bosonization 
and mapping to the quantum sine-Gordon model. 

In the quantum sine-Gordon model, neither operator 
d^cj) nor cos 6* can change the number of solitons. Thus, 
they can have non-vanishing matrix elements only be- 
tween the ground state and the states having an equal 
number of solitons and antisolitons and any number of 
breather. As a result, the Raman intensities predicted 
from the sinc-Gordon model in Sec. Ill Cl are in partial dis- 
agreement with the ones obtained from ED if the state S 
on Fig. |4] is identified with a sine-Gordon soliton. There 
are three different possibilities to resolve this discrepancy. 
The first one is that the state S is actually the second 
breather B2 of the quantum sine Gordon model. This 
is compatible with its symmetry since there is one such 
state in the -|-) sector. This is also compatible with 
the observed decrease of the weight of the S peak as D is 
decreased, as such a behavior is predicted by Eq. (|26|) . In 
that scenario, the ratio of the masses do not agree with 
the quantum sine Gordon prediction. The latter result, 
however, may result from having relatively large gaps 
that make the continuum limit unjustified. However, in 
this scenario, we have to understand why there is a state 
degenerate with the second breather in the sector (—,+). 

The second possibility is that the expression (PO)) . 
which is valid for Ja_/J\\ — » 00 has corrections of or- 
der {D/Jj_y"' with m an integer coming from higher 
order virtual processes. If these higher order correc- 
tion terms contain an operator such as (— 1)"'t^ or 
T^'^ the bosonized form of which are repectively cos 20, 



IV. CONCLUSION 

Although the theoretical interpretation of the Raman 
intensities obtained numerically at intermediate coupling 
in terms of the solitons and breather of the weak coupling 
sine-Gordon theory is still a partially open problem, from 
the point of view of experiment, we can safely conclude 
that in a two-leg ladder with Dzialoshinskii-Moriya in- 
teraction, the presence of a gap in the magnetized phase 
gives rise to bound states observable in magnetic Raman 
scattering. In the case of CuHpCl, earlier studies showed 
that typical parameters like J|| ~ 0.2 (as used in the nu- 
merics) and D ~ 0.05 (all in units of Jj_) give a fair de- 
scription of the magnetization curve^^^ and of the spe- 
cific heat in an applied magnetic field^^. Hence, the above 
investigation of Raman magnetic scattering (based on an- 
alytic and numerical techniques) clearly establishes that, 
for light polarization along the legs, at least two bound 
states should be visible, while at least a single bound 
state should be seen for light polarized along the rungs. 
We also expect that the spectral weight of the continuum 
predominantly appears at energies larger than typically 
^ 2J± i.e., in the case of CuHpCl, above ^ 3 meV (al- 
though the "theoretical" onset of the continuum is lo- 
cated just above the bound states as predicted by the 
SG description). In other words, a "pseudo-gap" should 
form between the bound states and ~ 2Jj_. In the par- 
ticular case of light polarized along the rungs, we further 
note that the bound state(s) alone should exhaust most 
of the overall spectral weight. We therefore believe that 
Raman scattering experiments should provide clear fin- 
gerprints of the existence of finite Dzialoshinskii-Moriya 
interaction in low-dimensional magnets. 
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